Thermofield Dynamics of Time-Dependent Boson and Fermion Systems 
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We formulate the thermofield dynamics for time-dependent systems by combining the Liouville- 
von Neumann equation, its invariant operators, and the basic notions of thermofield dynamics. The 
new formulation is applied to time-dependent bosons and fermions by using the time-dependent 
annihilation and creation operators that satisfy the Liouville-von Neumann equation. It is shown 
that the thermal state is the time- and temperature-dependent vacuum state and a general formula 
is derived to calculate the thermal expectation value of operators. 
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ly-j , The thermofield dynamics (TFD) introduced by Takahashi and Umezawa three decades ago is a canonical formalism 
for finite temperature theory to describe quantum systems in thermal equilibrium 0, • The great merit of TFD 
is that it preserves many useful properties of the zero-temperature field theory. The central concept in TFD is the 
thermal state, a pure state, in some extended Hilbert space, which corresponds to the thermal equilibrium, a mixed 
state, in the original Hilbert space. The three essential ingredients of TFD are (i) the tilde conjugation rule, (ii) 
the Heisenberg equation and (Hi) the thermal state conditions Q. In a simple physical terminology, the TFD of 
a quantum system doubles the degrees of freedom by introducing a fictitious Hamiltonian without any interaction 
Qh with the system and uses an extended Hilbert space of the direct product of the Hilbert spaces of the system plus 
the fictitious system. In the oscillator representation there is a temperature-dependent Bogoliubov transformation 
between the annihilation and creation operators of the total system and those of temperature-dependent ones. Then 
the thermal state is a two-mode squeezed vacuum state (temperature-dependent vacuum state) of the extended Hilbert 
space, which in turn is annihilated by the temperature-dependent annihilation operators. There have been numerous 
diverse applications of TFD to systems in condensed matter, nuclear physics, particle physics, quantum optics and 
cosmology in thermal equilibrium (for review and references, see Q and jj). 

On the other hand, there are also many systems exhibiting nonequilibrium characteristics. An open system in- 
teracting with a reservoir or a time-dependent system is such a nonequilibrium system. As the time-translational 
invariance of such a time-dependent system is broken, there is an obstacle in applying the Matsubara's imaginary-time 
method [j| . Neither does the conventional wisdom work well using the basis of the time-dependent energy eigenstates 
in evaluating thermal quantities since the density operator is not given by e~ l3H ( t K However, the closed-time path 
method by Schwinger and Keldysh is well-known and is widely used for such nonequilibrium systems 6J. Another 
canonical theory based on the (functional) Schrodinger equation has been employed to study nonequilibrium evolution 
of time-dependent systems |7[. Keeping the three ingredients of TFD, requiring the Hamiltonian to be tildian, and 
using the Heisenberg picture, Umezawa et al attempted to extend the TFD to such nonequilibrium systems [g. 

In this paper we show that when (ii) the Heisenberg equation, one of the ingredients of TFD, is replaced by the 
Liouville-von Neumann (LvN) equation, the TFD has a direct generalization to time-dependent boson and fermion 
systems, in particular, time-dependent oscillators whose mass and frequency may change explicitly in time. The 
invariant operators satisfying the LvN equation not only lead to the Hilbert (Fock) space of exact quantum states 
but also provide the correct density operator [T^j. Thus we are able to extend the TFD to time-dependent 
bosons (fermions) first by using the time-dependent annihilation and creation operators for the bosons (fermions) 
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and the fictitious bosons (fermions), all linear invariant operators [Tol ITU IT2I IT^. and then by finding the time- and 
temperature-dependent annihilation and creation operators through temperature-dependent Bogoliubov transforma- 
tion of TFD. Then the thermal state is a two-mode squeezed state of the time-dependent vacuum state for the bosons 
(fermions) plus the fictitious bosons (fermions). We find a general formula for evaluating the thermal expectation 
values of operators and finally discuss the distribution which evolves from an initial boson distribution. 

The organization of this paper is as follows: In Sec. II, we briefly review the TFD for boson and fcrmion systems. 
In Sec. Ill, we introduce the time-dependent annihilation and creation operators for a time-dependent boson. In Sec. 
IV, we extend the TFD to time-dependent boson and discuss the physical implication of the TFD for a boson which 
evolves from an initial thermal state to a final one through a time-dependent interaction. In Sec. V, TFD is extended 
to time-dependent fermions. Additional comments and conclusion are given in Sec. VI. 

II. TFD FOR BOSON AND FERMION SYSTEMS IN EQUILIBRIUM 

We briefly review the TFD for static (time-independent) bosons and fermions in a way that can be readily applied 
to time-dependent ones in the following sections. The static boson (fermion) has the Hamiltonian of the form 

H = fuva^a, (1) 
where the standard annihilation and creation operators satisfy the commutator (anticommutator) 

[<z,at]± = l. (2) 
Quantum statistics of the boson (fermion) in thermal equilibrium is described by the density operator 

P = ^ H , Z = Tr[e-? H ], (3) 

where (3 is the inverse temperature. For instance, the thermal expectation value of an operator A is given by 

(A) T = Tr[pA] = i £ ^ n {n\A\n), (4) 

n 

with \n) being number states. The thermal equilibrium is a mixed state with the probability p n — e~P hn for each 
projector \n)(n\. 

The idea of TFD is to double the system by adding a fictitious system and extend the thermal equilibrium in the 
system's Hilbert space to a thermal state, a pure state, in the extended Hilbert space of the total system. For that 
purpose, let us introduce a fictitious boson (fermion) using the tilde conjugation rule, with the Hamiltonian, 

H = fttoa^a, (5) 

with the number state \h) and the commutator (anticommutator) 

[a,a)]± = l. (6) 

As {a, a^} and {a,a^} describe two independent systems, H and H, respectively, they commute (anticommute) with 
each other 

[a,a]± = [a, a f ]± = [a T ,a]± = [o^a^i = 0. (7) 

The total Hamiltonian 

H := H - H = Huj(a'a- a f 5) (8) 

now carries an extended Hilbert space H = H <g> H. Each state of the extended Hilbert space consists of a product of 
number states 

\n, rh) = \n) ® |m). (9) 
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In TFD the thermal equilibrium corresponds to the thermal state in the extended Hilbert space, described by a 
pure state 



|O09)> = 7^E e "^ n/2 i at ' l(it ' l l ' 5 ) 



ZV2 — n \ 

n 

= ( 1=Fe -/5")l/2 e -(/3^/2)atat| ^ ^ 

where the upper (lower) sign is for bosons (fermions) and |0) = |0, 0) in the second line. In fact, the thermal state is 
a two-mode squeezed vacuum state obtained by applying a Bogoliubov transformation to the vacuum state 

|0(/3)}=e- G |0), (11) 

where 

G=-i6{P){aa-a)at). (12) 
Here (3(9) is a temperature-dependent parameter determined by 

cosh#(/3) = (l-e^ 1 ")- 1 / 2 , 

sinh<?(/?) = e -^/2(! _ e -/s««)-i/a j ( 13 ) 

for bosons and 

cosfl(/3) = (l + e-^)- 1 / 2 , 

sin0(/3) = e -/3W2( 1 + e -/3^)-i/2 ) ( 14 ) 

for fermions. 

The two-mode squeeze operator introduces the temperature-dependent annihilation and creation operators through 
the Bogoliubov transformations 

a(f3) = eT lG ae lG = cosh 6>(/3)a - sinhflO^a 1 ', 

a(/3) = e _tG ae lG = cosh0(/3)a - sinhtfCS)^, (15) 

for bosons and 

o(/3) = e- lG ae lG = cos6(/3)a-sin9(0)a\ 

a(j3) = e~ lG ~ae lG = cos (9(^)5 + sin ^(^a 1 , (16) 

for fermions. Similar relations follow for o' (/3) and a'(f3) for bosons and fermions. The thermal state in Eq. Ijllfl is 
nothing but a temperature-dependent vacuum state 

a09)|O(/3)> = o(/9)|0(/3)> = 0. (17) 

Now the expectation value of an operator A with respect to thermal equilibrium is equivalent to that with respect to 
the thermal state 

(A)i = (0(/3)|A|0(/3)). (18) 
The expectation value of the number operator leads to the boson distribution 

(0(/3)|at a |0(/3)) = sinh 2 6(/3) = ^ _ v (19) 



and the fermion distribution 



(0C9)|oto|009)) - sin 2 9{P) = ^ +v (20) 
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III. TIME-DEPENDENT BOSON SYSTEM 



We now consider a time-dependent boson system and find the ingredients necessary for TFD extension. The most 
general time-dependent quadratic Hamiltonian for bosons takes the form 



1 



1 



H(t) = h oo (t)a t a + -u) + (t)a u + -cola 



(21) 



where a and are the Schrodinger (time-independent) annihilation and creation operators, and uiq is real and lo + 
real or complex. In particular, the oscillator with time-dependent mass and frequency, 



H(t) = 



m(t) 



2m(t) 



(22) 



belongs to the Hamiltonian, Eq. 121|) . Note that p and q in Eq. 11221) are also the Schrodinger operators, but 
the Hamiltonian depends explicitly on time through m and uj. An instantaneous energy eigenstate \e n ,t) of the 
Hamiltonian, Eq. I|21l) . is definitely not an exact quantum state of the Schrodinger equation. Moreover, the density 
operator is not given by e~$ H ^) . So the set {|e„,t)} is not a good basis to study the noncquilibrium evolution of the 
thermal system. 

We explain a physical motivation for replacing the Heisenberg equation in the ingredient (ii) by the LvN equation. 
In the Heisenberg picture, it is the Heisenberg equation that determines the operators evolving in time. In fact, in 
terms of the evolution operator U of the Schrodinger equation, each Schrodinger operator A leads to the Heisenberg 
operator Ah — W AU . For a time-dependent system, these operators carry all the quantum information of the system. 
As states do not evolve in time, the initial equilibrium with the density operator p leads to the thermal expectation 
value of A given by 



(A) T = Tr[pA H (t)] = Tt[UpU^A}. 



(23) 



Therefore, to know thermal properties of the system, either the Heisenberg operator Ah or the density operator 
p(t) = UpW should be given in advance. 

On the other hand, there exist the invariant operators that satisfy the LvN equation [9| ■ A complete set of invariant 
operators provide another picture for a time-dependent system and such a set of invariant operators are known for 
time-dependent oscillators in terms of classical solutions 0, 0] . In this sense the LvN equation may replace the 
Heisenberg equation in TFD. Further, the linearity of the LvN equation allows the construction of the density operator 
from an invariant operator |ldj |. In fact, we may find the time-dependent annihilation operator, an invariant operator, 
of the form 



a(t) = /(-)(t)o + /C+)(t)ot 
and its Hermitian conjugate a'(t), and impose the LvN equations 

ih^ + [a(t),H(t)]-=0, 
at 



(24) 



(25) 



The time-dependent creation operator a' (t) is another invariant operator. The pair {a(t), a' (t)} form a complete set. 
Therefore any invariant operator can be constructed out of them. 
The LvN equation (|25() leads to the vector equation 



where 



and 



dV 

i— + M V V = 0, 
at 



(26) 



(27) 



M V {t) = CJ (7 3 - ~{UJ* + - W+)<71 - ~{UJ* + + LJ+)<J 2 - 



(28) 
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Here a's are the Pauli spin matrices. The solution to Eq. (|26|l provides the time-dependent annihilation and creation 
operators. Further, the equal-time commutator 

[o(t),ot(t)]_ = l (29) 
can hold by choosing the initial data V^azV = 1 since Eq. (|26[1 leads to the relation 

|(FV 3 F)=0. (30) 

For the time-dependent oscillator, Eq. (|22J) . the following form of the time-dependent annihilation operator is known 
dill 

a(t) = ^=[v*(t)p-m(t)v*q]. (31) 

The time-dependent creation operator a*(t) is the Hermitian conjugate of a(t). In fact, these operators are invariant 
operators satisfying the LvN equation (|25|l for the Hamiltonian, Eq. I|22|l . when v is a complex solution to the classical 
equation of motion 

m + ^rrMt) + u 2 (t)v(t) = o, (32) 

m(t) 

and satisfies the Wronskian condition 

m(v*v — vv*) = i. (33) 

The number states of the number operator N(t) — a^(t)a(t), another invariant operator, defined as, 

N(t)\n,t) =n\n,t), (34) 

are exact quantum states of the Schrodinger equation and constitute the Hilbert space 0, 0, 0, . The state 
|0,i) is the time-dependent vacuum that is annihilated by a(t), and the number states are obtained by applying the 
creation operators on it 

\n,t) = —L>\0,t). (35) 
Vn! 

It follows that the density operator may be given by 

Z 

where (3 and oj are constants. In the case of the oscillator, Eq. 112211 . the position and momentum operators have the 
oscillator representation 

q = Vh[v(t)a(t)+v*(t)a^(t)], 

p = Vhm(t)[v(t)a(t) +v*(t)a?(t)]. (37) 

The merit of using the LvN equation is that we can keep all the steps of the time-independent boson case in finding 
the Hilbert space, the density operator and so on. 

IV. TFD FOR TIME-DEPENDENT BOSON SYSTEM 

To extend TFD to the time-dependent boson system, Eq. JHJ, we use the tilde conjugation rule (cA) = c*A and 
introduce a fictitious boson system with the Hamiltonian 

H(t) = h uj Q {t)a*a + ^uj* + {t)5,V + ^+a 2 . (38) 
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As for the time-dependent boson case, we introduce the time-dependent annihilation operator for the fictitious boson 

a(t) = /<->*(t)o + /<+>*(t)ot ) (39) 
where f^' satisfy Eq. Ij26[l . Then a(t) and its Hermitian conjugate a^(t) satisfy the LvN equations 

iH^ + [a(t),-H(t)}-=0, 
ift^£ + [a t (t),-.ff(t)]-=0. (40) 

The equal-time commutator also holds 

[5(*),at(t)]_ = i. (4i) 

The number states 

N(t)\n, t) = a f (t)a(t)\h, t) = n\n, t) (42) 

are the exact quantum states for the Hamiltonian, Eq. (|38ll . 
The total Hamiltonian is now given by 

H(t) = H(t)-H(t), (43) 
where the operators of the boson and their tilde operators commute with each other 

[a, a]_ = [a, a f ]_ = [a f , a]_ = [a f , a]_ = 0, (44) 

and 

[a(t), a(t)]- = [a(t), t(t)]_ - [tf(t), a(f)]_ = [a+(0, a(t)]- - 0. (45) 

In fact, a(i), a'(t), a(i) and a^(t) are the invariant operators satisfying the LvN equations for the total Hamiltonian, 
Eq. (|4*3l . The Hilbert space of the total system consists of 

|n,m,t) = |n,t) <g> |m,t) = J_ ; |0,0,t). (46) 

v«! vm! 

The density operators may be defined by 

p(t) = ie-^IW), (47) 

z 

p(t) = i e +^ st W a W, (48) 

which obviously satisfy the LvN equations. Here /3 and w are constants that may be fixed by the initial temperature 
and frequency. The density operator in the extended Hilbert space is given by 

Pit) = pit) ® pit) = J- e -Pr*o("Ht)a.(t)-aHt)a(t)) (49) 

The thermal expectation value of the operator A of the system now takes the form 

(A) = Tr[p(t)A] = (0(/3), t|A|0(/?), t), (50) 
where the thermal vacuum state is given by 

W),t) = _L_^e-^"/ 2 latn ( ^ W |o,o,t) 

n 

= y/l - e -/5^ e -(/3^/2)at( t )a+(t)| 0i ^ 5 (51 ) 
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with |0, t) = |0, 0, t). The thermal state is an exact eigenstate of the Schrodinger equation for the total system, Eq. 
(|43|l . The thermal state is also written as a time-dependent two- mode squeezed state of the vacuum state 

\O(0),t) =e- iG W|0,t>, (52) 

where 

G(t) = -i6(0)[a(t)a(t) - at (t)tf (t)]. (53) 

Here 6(f3) is the same parameter fixed by Eq. I|13|) . 

As the density operator in Eq. (|36|l or (|48|l involves a constant 0, we may find the time- and temperature-dependent 
annihilation and creation operators through the Bogoliubov transformation 

a(0,i) = cosh9(0)a{t) - smh6(0)a! (t), 

a(0,i) = cosh9(0)a(t)-sinh9(0)a^(t), (54) 

and their inverse transformation 

a(t) = cosh9(0)a(0,t) + smh8(P)a!(0,t), 

a(t) = cosh9(0)a(0,t) + sinh6»(/3)a t (/3,t). (55) 

We get similar equations for a'(0,t), cJ{(3,t), a'(t) and 5^(t) by using the Hermitian conjugate of these equations. 
As 9(0) is a constant, a(0,i), a(0,t) are invariant operators. Then the thermal state is the time- and temperature- 
dependent vacuum 

a(0,t)\O(0),t) = a(0,t)\O(0),t) = O. (56) 

The thermal state \O(0),t), as an eigenstate of the invariant operators a(0,t) and 5(0, t), is an exact eigenstate of the 
total system. At each moment, the boson still keeps the same boson distribution since the expectation value of the 
time-dependent number operator yields 

(0(/5),t|at(t)a(f)|0(/?),i) = sinh 2 ^(/3) - ^ _ v (57) 

Using TFD for time-dependent bosons, we are able to find the thermal expectation values of operators. In general, 
through the Bogoliubov transformations from {a(t), a'(t)} to {a(0, t), a>(0, t)}, we find the formula 

(F(a(t),a f (t))) T = (0(0), t\F(cosh9(0)a(0, t) + sinh 9(0)5^ (0, t), 

cosh6(0)(J(0,t) + smb. 9(0)5,(0, t))\Q(J3),t}. (58) 

This provides the basic rule for calculating matrix element of any operator in TFD. For instance, in the case of the 
oscillator, Eq. Ij22(l . using the position representation in Eq. (|37[1 

q = Vhcosh9(0)[v(t)a(0,t) +v*(t)a)(0,t)} 

+Vhsinh9(0)[v*(t)a(0,t) + v(t)a t (0,t)\, (59) 



we obtain 



(O(0),t\q 2n \O(0),t) = H n J2 (lf}(O,t\cosh 2k 9(0MtM0,t)+v*(t)aH0,t)] 2k 



k=0 

x sinh 2 "- 2 * 8(0)[v*(t)~a(0, t) + v(t)a*(0, i)] 2 "~ 2fc |0, t). (60) 
After normal ordering, we finally obtain the result 

(q 2n ) T = (O(0),t\q 2n \O(0),t) 



(2n)\ 
2 n n\ 



[Hv* (t)v(t)] n (1 + 2 sinh 2 9(0)) n . (61) 



We discuss the physical implication of TFD for a time-dependent boson, when it evolves with a time-dependent 
interaction from initial lo^s at t = ti to final ones, w/'s at tf. That is, all w's change from Wj's to w/'s. The solution to 
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Eq. (|26|l is necessary in finding the time-dependent annihilation operators in Eqs. I|24|) and l|39|l and their Hermitian 
conjugates. However, from constants wy's, we may find a Bogoliubov transformation of the form 

ai = ficif + va}j, 

a\ = ii*a^ + v*af, (62) 

where {ai,a\} for uj^s and {aj,a^} for w/'s. Here n and v, which should be determined by the solution to Eq. 
(|26|l , carry all the information about the history of interaction and may take the form 

fi = fi(ti,tf,ci)i,Wf), v = v(ti,tf,uji,u>f), (63) 

and satisfy 

(i* fi — v*v = 1. (64) 

If the boson is initially in thermal equilibrium with the inverse temperature (3 and has the boson distribution hi = 
l/(e^ Sw — 1), then, according to Sees. Ill and IV, the boson in the final state has a different distribution 

(0(P),t f \al ai \0(P),t f ) = v*v+ \^2 V [\ - (65) 

The first term is originated from the particle production from vacuum fluctuations |l4j . (0, tf |a|aj|0, tf) — v*v, and 
the second term is a purely thermal result, having an overall amplification factor (1 + 2v*v) to the boson distribution. 
Thus the evolution of the time-dependent system leads to a distribution quite different from the boson distribution 
function. 



V. TFD FOR TIME-DEPENDENT FERMION SYSTEM 



The time-dependent fermion system, quadratic in the annihilation and creation operators, has the Hamiltonian 

H(t) = ft[wo(0(a f a - tfb) + u + (t)aW - uo* + {t)ab + uj_{t)atf - w* (t)a% (66) 

where a, for the particles and b, b^ for the antiparticles are all the Schrodingcr (time-independent) operators, and 
wo is real and w± are real or complex. They satisfy the anticommutators 



a, a 



[b, = 1, [a, b} + = [a, &t] + = [ a t, &t] + = . (67) 



The Hamiltonian, Eq. H66|l . is a Hermitian operator and its time-dependency comes only from the parameters u>'s. 
As for the boson case, we may find a pair of time-dependent invariant annihilation operators |T^ 

a(t) = /(-)( i )a + /(+)(t)at+ 5 (-)(t)6 + 5 (+)(t)6t ) 

b(t) = ft\t)a + A + \tW +g[- ) {t)b + g[ + \t)b\ (68) 

and the invariant creation operators a) (t) and w (t) are the Hermitian conjugates of a(t) and b(t), respectively. The 
LvN equations for the operators in Eq. l(6"5|l with the Hamiltonian, Eq. (|Bl)jl . lead to the following vector equations 

m 

dW 

i^— + cuoCTiW + M W Z = 0, 
at 

i— + u aiZ + M Z W = 0, (69) 

where 

1 ( fH + f (+ 
W(t) = -= 1 1 +J 



1_ (g(-)+gM 
V2 \ 9 { - ] ~ 9 { 
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and 



1 1 i 1 

M w{t) = _ -(u>*_ +uJ-)ai - -(w* +cj+)o- 2 - -^+) cr 3, 

1 . 1 .. , i, . . 1 



M z {t) = --(w* +w_) ( 7 1 + -( w ; + w + )(7 2 + -(w;-w + ) ( T3, (7i) 

where / is the identity matrix. The equal-time anticommutators 

[o(t), ot(i)]+ = W),b\t)] + = 1, [ (i), &(<)]+ = [a(f), 6t(<)] + = [ a t(i), fct( t )] + = o, (72) 
are guaranteed by the relations 

^(W^W + Z^Z) = 0, 

J^> 3 ^ + Zja 3 Z b ) = 0, 

f^ttf^ + ^Z b ) = 0. (73) 

Each set of {W a , Z a } and {Wb, Zb] gives rise to a(t) and 6(i), and their Hcrmitian conjugates. The number operators 
N a (t) — a'(t)a(t) and Nb(t) — b'(t)b(t), which are also invariant operators, span the state vector space of the fermions 

\0,t), ot(t)|0,t), tf(t)\0,t), ot(t)6t(t)|0,t). (74) 

Here the vacuum state, |0,t) = |0, 0, t), implies no particles and antiparticles. The density operator may take the 
form 

p M = jL e -^(otWo(t)-6t(t)6(*)) ) (75) 
Z 

where f3 and ui are again constants. 

To construct the TFD for fermions, we introduce the fictitious fermion Hamiltonian 

H(t) = h[uj a (t)(a^a - ftb) + w$.(t)aty - uj + (t)ab + u*_(t)aP - w_(i)a t 6]. (76) 

Then there are the time-dependent annihilation operators, invariant operators, 

a{t) = fi->{t)~a + fi + >{t)~a^+g^>{t)b + g^>{t)b\ 

b[t) = / b ( - ) *W5 + / b (+) *(t)at +ff (-)*(t)6 + ff (+)*(t)6t. (77) 

These operators satisfy the LvN equations (|40|l when f's and g's satisfy the vector equations H69|l . We are then 
equipped with the time-dependent annihilation and creation operators for the total system 

H(t) = H(t) - H(t). (78) 

The time-dependent vacuum state of the total system is annihilated by all the annihilation operators of the fermions 
and the fictitious fermions: 

a(t)\0,0,t) = b(t)\0,0,t) = a(t)|0,0,i) = b(t)\0,6,t) = 0. (79) 

Here the vacuum state, |0,0, t) — |0, 0,0,0,4), implies no particles or antiparticles and their counterparts. However, 
we will continue to use the abbreviation |0, 0, t) for this state. 

The thermal state of TFD is the two-mode squeezed state of the time-dependent vacuum state 

\Q(fi),t) =e- lG ^|0,i), (80) 

where \0,t) = |0, 0, t>, and 

G F (t) = -i6(j3)(aa - at at +bb- b r tf). (81) 
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Here 6(P) is the same parameter fixed by Eq. (|14(l . The two-mode squeeze operator introduces the time- and 
temperature-dependent annihilation and creation operators for the particles through the Bogoliubov transformation 

a(/3,t) = cos9{P)a(t)-sin9(P)a^(t), 

a(/3,t) = cos9(p)a(t) + sin 9 (P)a f (t), (82) 

and for the antiparticles 

b(/3,t) = cos9(p)b(t)-sin9(p)tf(t), 

b(P,t) = cos 0(/3)b{t) + sm6(p)tf{t). (83) 

We get similar equations for a^(P, t), a^(0, t), t) and w((3,t) by using Hermitian conjugate on these equations. 
Their inverse transformations are 

a(t) = cos9(P)a{P,t) +sin0(P)a}(P,t), 

a(t) = cos9(P)a(P,t)-sin0(P)al(P,t) (84) 

and 

b{t) = cos 9((3)b(f3,t) +sm9(P)P(P,t), 

b(t) = cos9(P)b(P,t)-sin9(p)tf(p,t). (85) 
Further, the thermal state is annihilated by 

a(p,t)\0(p),t) = a(p,t)\Q(p),t) = b(p,t)\0(p),t) =b(/3,t)\0(j3),t). (86) 
Now the expectation value of the particle number operator 

(0(p),t\a\t)a(t)\0(p),t) = sin 2 9=—^. (87) 

The thermal expectation value of operators, for instance, of the particles is given by the formula 

(F(a(t),a*(t))) T = (0(f3),t\F(cos9((3)a(/3,t) +sin0(/?)a t (/3,t), 

cos6»(/3)a t (/3,t) + sin0(P)a(P,t))\O(P),t). (88) 



VI. CONCLUSION 



To summarize, in this paper we have used the time-dependent annihilation and creation operators of the LvN equa- 
tion to complete TFD for time-dependent boson and fermion systems. The first ingredient, (i) the tilde conjugation 
rule, has been accomplished by introducing the fictitious boson or fermion operators and appropriately constructing 
the extended Hilbert space. The second ingredient, (ii) the Heisenberg equation, is replaced by the LvN equation. 
The annihilation and creation operators from the LvN equation have led to exactly the same procedures as for TFD 
of time-independent boson and fermion systems. The third ingredient, (iii) the thermal state condition, is guaranteed 
by the time- and temperature-dependent Bogoliubov transformations, Eqs. 1|54|1. (|82|l and l|83() . The thermal state 
for time-dependent bosons is the time- and temperature-dependent vacuum as in Eq. I|56|l . which can be written as 

a(t)\0(P),t) =tanh6'(/3)a t (t)|0(/3) ! t), a(t)\0(/3),t) = tanh 6(p)<J(t) |0(/3),t). (89) 

The thermal state condition for time-dependent fermions is given by 

a(t)\0(p),t) = tan9(P)tf(t)\0(P),t), a(t)\0(P),t) = -ttm9(P)at (t)\0(P),t), 

b(t)\0(p),t) = tan9(Pp(t)\0(p),t), b(t)\0(p),t) = -tan 9(p)tf(t) \0(P),t). (90) 

Note that the thermal state conditions, Eqs. I|89|) and (|90|l . still have the same form as for the time- independent 
boson and fermion systems. We may thus conclude that the LvN equation provides a direct generalization of TFD 
to nonequilibrium systems such as time-dependent boson and fermion systems. Finally it should be stressed that the 
replacement of the Heisenberg equation by the LvN equation leads us from the Heisenberg picture to the Schrodinger 
picture. Such a procedure yields interesting and useful results for quantum mechanical time-dependent systems. 
Extension to quantum field theory may be possible and this would bring the open systems to a treatment similar to 
the case of systems in equilibrium. That would provide a true extension of TFD, a la Takahashi and Umezawa that 
has proved so very successful for treating equilibrium problems, to systems out of equilibrium. Hopefully this will 
provide a useful perspective on such a class of problems. This topic is under active consideration. 
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